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COMPARISON OF ONE- AND TWO-DIMENSIONAL HEAT-TRANSFER 
CALCULATIONS IN CENTRAL FIN-TUBE RADIATORS 
byNorbertO. Stockman, Edward C. Bittner, and Earl L. Sprague 

Lewis Research Center 

SUMMARY 

An analysis is given of the two-dimensional heat transfer, including gray body radiant 
interchange, in the cross section of a central fin-tube radiator panel. Results of this 
analysis are used to evaluate several one- dimensional methods of varying complexity for 
calculating the heat rejection rate of a central fin- tube radiator panel. Most methods 
gave good agreement with the two-dimensional results. In view of the excellent agree- 
ment of one of the simpler methods which neglects tube wall temperature drop and ac- 
counts for radiant interchange between fin and tube simply by using the projected area of 
the tube, it seems unwarranted to use the more complex methods which gave no better 
agreement. Details of the numerical method of solution of the two-dimensional equations 
are given in an appendix. 


INTRODUCTION 

Electric power generation systems for space applications must reject large amounts 
of waste heat. At the present state of the art, the most likely method of rejecting this 
heat is by means of a radiator that utilizes some sort of fin and tube configuration. Var- 
ious fin and tube combinations have been proposed, but only the central fin-tube geometry 
(fig. 1) will be considered in this report. Since the radiator represents a large portion 
of the powerplant weight, it must be accurately designed. Radiator design calculations 
are usually based on the assumption of one-dimensional heat transfer in the cross section 
of the fin and tube panel. Several different methods have been used to compute the heat 
rejection from central fin- tube panels (e. g. , refs. 1 to 8). These methods contain vary- 
ing degrees of complexity, and it is not obvious which methods are more accurate. 

In order to evaluate the one- dimensional methods, two-dimensional calculations of 
fin-tube heat transfer were carried out at NASA Lewis Research Center. These calcula- 
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Figure 1. - Fin-tube radiator. 


tions involved fewer assumptions than any of the one- dimensional methods and therefore 
are assumed to give a more accurate representation of the heat rejection. The analysis 
of the two-dimensional conduction in the cross section of a central fin-tube radiator panel 
with radiation from both surfaces and radiant interchange among the fin and tube surfaces 
is presented herein. The numerical method of solution is detailed in an appendix. Heat 
rejection obtained by the two-dimensional calculations is compared to that obtained by 
several different one-dimensional methods for 39 specific design geometries divided into 
the following three groups: (1) 13 miscellaneous designs for various power system appli- 
cations; (2) 11 cases arbitrarily created to study the effect of changes in the fin- tube de- 
sign variables on the comparison between one- and two-dimensional results; and (3) 15 off- 
design configurations based on a near-minimum-weight design for a Brayton cycle radi- 
ator. 

ANALYSIS 

The heat rejection rate of the cross section of the fin-tube panel of a radiator of the 
type shown in figure 1 will be computed by several different methods. Details of the two- 
dimensional analysis, including the calculation of the heat rejection rate, will be given, 
and formulas for calculating the heat rejection rate by several one-dimensional methods 
will be presented without detailed development. 
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Assumptions 


All the methods (one- and two-dimensional) are based on steady- state heat conduction 
in the cross section of the fin and tube, radiation from both fin and tube surfaces, and the 
following specific assumptions (other more restrictive assumptions are required for cal- 
culating heat rejection rate by the various one-dimensional methods, and these will be 
given later with the formulas for heat rejection rate for each method): 

(1) The radiator is infinitely long and made up of an infinite number of identical finned 
tubes. 

(2) The inside tube wall is isothermal both longitudinally and circumferentially. 

(3) The tubes and fins are made of the same material, and the material properties 
are assumed invariant with temperature. 

(4) If there is a liner (as in fig. 1), it is of the same material as the tube and is in 
perfect thermal contact with the tube. 

(5) Incident radiation from external sources such as Sun and planets and adjacent ve- 
hicle components is accounted for by a completely encompassing surface of constant tem- 
perature called the equivalent sink temperature T„ (ref. 9). This temperature is as- 
sumed to be the same on both sides of the radiator. 

(6) Radiating surfaces are gray at a prescribed value of emittance invariant with tem- 
perature. 

(7) Absorptance is equal to emittance and equal to one minus the reflectance. 

The assumptions of the same sink temperature on both sides of the radiator and the 
same inside temperature for all tubes assures symmetry about the three center lines of 
the section as shown in figure 2, and thus only the shaded portion need be analyzed. Radi- 
ant interchange with surfaces outside the shaded region is included, however. It should 
also be pointed out that, in the two-dimensional analysis, the outside tube surface is not 
isothermal, and radiant interchange between individual points (increments of area) on the 
fin and tube surfaces is taken into account. Furthermore, this interchange extends into 
the axial direction (perpendicular to the plane of the cross section) and is accounted for in 
the development of the interchange view factor (see next section). 

Since the one-dimensional methods require all the assumptions of the two-dimensional 
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Figure 2. - Two-dimensional cross section of fin-tube radiator showing region analyzed (shaded area). 
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plus additional ones, the two-dimensional results are the limit to which the one- 
dimensional results tend with increasing refinement. Thus the two-dimensional heat re- 
jection is used as the standard to evaluate the relative merits of the one-dimensional 
methods. In the next section, the differential equation and boundary conditions governing 
the two-dimensional heat transfer in the fin cross section under the assumptions listed 
previously will be presented. 


Equations and Boundary Conditions 

The steady- state two-dimensional heat flow in the fin and tube cross section is given 
by 


V^T = 0 (1) 

where T is the temperature. (All symbols are defined in appendix A. ) The Laplacian 
2 

operator V is written in polar coordinates in the tube and in rectai^lar coordinates 
in the fin (fig. 3) to simplify grid generation for the numerical solution. 

Thus, in the tube 


1-5I + 

r 0r 



J_ ^ 

r^ dtp^ 


= 0 


and in the fin 

d^T ^ d^T ^ Q 

.2 
ax dy 


y 



( 2 ) 


( 3 ) 
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These are the two partial differential equations which must be solved. The solutions must 
also satisfy the following boundary conditions. Along the inside tube wall at r = r^^^, the 
temperature is prescribed; 


(4) 

There is no heat flow across the lines <p = 0 and ^ = (ir/2) in the tube because of sym- 
metry, so 

(il\ .0 (5) 

Also, because of symmetry, there is no heat flow across the lines y = 0 and x = f + r^ 
in the fin, so that 



On the surface of the tube at r = r^ there is radiation to space where the sink tem- 
perature is T- and interchange with surfaces of the fin and the adjacent tube. The 

s 

boundary condition for net heat rejection from the tube surface can be written as 

-k/^ ) = eaT"^ - q!H = € (aT^ - H) (7a) 

where e is the emittance, a is the absorptance of the surface, and H is the incident 
radiation including that from the sink and from other parts of the raditor surface. Simi- 
larly on the surfe.ce of the fin at y = t, the net heat rejection is 

_k = e (oT'^ - H) (7b) 

9y 

Since that part of the incident radiation H contributed by the radiator surface is unknown 
and must be determined as part of the solution of the problem, H must be written in a 
form in which the radiator surface radiation appears explicitly. The equation for H, 
derived in appendix B, is 
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( 8 ) 



where is the radiosity (the total energy per unit time and unit area) leaving an ele- 

ment dA on the fin-tube surface and for a particular element dA* is given by 

®dA* “ ^ ^^dA* ^®dA* 

Substituting equation (8) into equation (9) and replacing p by 1 - e yield 



In equations (8) to (10) is the view factor from a particular element dA* 

to any other element dA; that is, it is that fraction of the total heat per unit time leaving 
dA* which strikes dA. The integrals in these equations represent integrals over all 
parts of the fin and tube visible from dA*. Thus, for example, if dA* is on the left 
tube (fig. 4), then 


/ /^U /*‘^2U 


rod<^l-rod<P2 


( 11 ) 


where dF^, factor from the element r^ d(p^ on tube 1 to the element 

dx on the fin, Xj^ and x^ are the lower and upper limits, respectively, of visibility on 
the fin from v^dicp^, dF^ is the view factor from element r„ d.<p^ on tube 1 


dF„ T. 1® the view factor from element r dw , 
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to element dcp 2 on tube 2, and (pgL *^2U lower and upper limits, re- 

spectively, of visibility on tube 2 from r^ d(p (View factors and limits of visibility for 
several fin-tube radiator configurations are given in ref. 10. ) 

It should be noted that the areas here (and heat rejection rates later) are per unit 
' axial length. This is to avoid infinite quantities since the axial length is assumed infinite. 
The view factors, however, are actually between elements of area of infinite length. (This 
presents no problem since view factors are fractional quantities and may be defined for 
infinite areas. ) For example, the view factor dF ^ is that portion of the radiation 

leaving the area r^ dq>.^ by infinite length on tube 1 and strikii^ the area dx by infinite 
length on the fin. It can be seen that the radiation leaving the surfe.ce at a certain axial 
location on the tube but striking a different axial location on the fin is accounted for in the 
view factor. 

Equations (2) to (8) and equation (10) in addition to the prescribed geometry r r^, 

L, and t and the temperatures T. and T„ determine the temperature distribution 
throughout the fin and tube. Equations (2) and (3) are differential equations in T; equa- 
tion (10) is an integral equation in B; and T and B are related through equations (7) 
and (8). The equations were solved numerically on an IBM 7094 using a finite difference 
block overrelaxation method. Details of the numerical method of solution are given in 
appendix C. 


Heat Rejection 


The temperature distribution resulting from the solution of the equations of the pre- 
ceding section is used to obtain a net heat rejection rate. The net heat rejection rate is 
also calculated by several one -dimensional methods for comparison with the two- 
dimensional results. 

Two-dimensional . - Once the two-dimensional temperature distribution has been ob- 
tained, the net heat rejected Q2JJ by the radiator per unit axial length for one quadrant 
can be obtained by integrating the left-hand side of equation (7a) over the tube outer sur- 
face and the left-hand side of equation (7b) over the fin outer surface; that is. 


Z <P=(7t/2) 

sin“^(t/rj 


Q. 


2D 



( 12 ) 


One-dimensional. - Several one-dimensional methods of calculating the net heat re- 
jection rate will be considered. The formula and a brief description of each method. 
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which will hereinafter be referred to by its number, follows. All methods are based on 
the assumptions of constant tube outer surface temperature and of one-dimensional heat 
conduction in a radiating fin. The conduction in the fin is accounted for by a fin efficiency 
or an overall efficiency. The fin efficiency depends, in general, on two dimensionless 
parameters: the conductance parameter X = coT^f /kt, and the sink temperature ratio 
Tg/T|j. The temperature T^^ is the fin base temperature and is equal to T^ or T.^^ 
depending on, whether or not the tube wall temperature drop is taken into account. When 
radiant interchange with the tube is taken into account, an overall efficiency is used that 
depends on X and Tg/T^^ and also on a third parameter, the fin-tube profile ratio f/r^. 

The items that differ from method to method are whether or not a radial one- 
dimensional temperature drop across the tube is considered (i. e. , whether Tb equals 
Tq or T.jj); how radiant interchange between the radiator surfaces is accounted for, if 
at all; what area is used for the tube radiating surfece; and how nonblackbody effects (e) 
are handled. In all methods except (6), the nonblackbody effect is handled simply by in- 
troducing € into the formula for Q. The methods, which follow, are summarized in 
table I. 

Method 1: 


(13) 

where is the fin efficiency based on a fin base temperature equal to T-^. Fin effi- 
ciency for this and the next three methods includes the effect of the sink temperature T_ 

s 

but not the effect of interchange with the adjacent tubes. 

No temperature drop in the tube is considered (i. e. , T^ = Tj^), radiant interchange 

is neglected, and the actual area of the tube is used. The fin efficiency is a function 

primarily of X. and secondarily of T_/T. and can be obtained from reference 9. 
in s in 

Method 2 : 


This is the same as method 1 except that the tube projected area is used instead of the 
actual area. 
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TABLE I. - SUMMARY OF ONE -DIMEJ^KIONAL METHODS OF CALCULATING 


HEAT REJECTION RATE Q 


Method 

Formula 

Temperature 
drop across 
tube wall 

Interchange 

Tube area 

1 


No 

None 

Actual 

2 


No 

None 

Projected 

3 

Qs = - Ts)(V + 

Yes 

None 

Projected 

4 

Q4 = 0.85 

Yes 

Approximated by 

Actual 


' ' ' is ' 


factor of 0. 85 


5 

Qs = - Ts)(^ + ^o) 

Yes 

Accounted for in 

Projected 




7]rj^ based on 


6 

% = vK’^0 - Ts)(^ + V 

Yes 

Accounted for in 

Projected 




?7* based on N^ 



Method 3: 


< 15 ) 

This is the same as method 2 except that a temperature drop across the tube wall is con- 
sidered and the fin efficiency tj is based on fin base temperature equal to the tube outside 

wall temperature and is a function of and T„/T^. 

o s' o 

Method 4: 

Q4=0.8fcc-(T^-T^)(V + 

This method uses the tube wall temperature drop and the actual area of the tube and ac- 
counts for the radiant interchange by the approximate interchange factor of 0. 85. The 
results of reference 8 are based on this method. 

Method 5: 



( 17 ) 
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This method uses the temperature drop in the tube and the projected area of the tube. 
Both the radiant interchange and the fin temperature distribution are accounted for by an 
overall effectiveness rj,p which is a function of X^, Tg/T^, and f/r^ (see, e. g. , refs. 4 
and 6). 

Method 6: 


Q6=»^<’(T^-T^)c«+ro) (18) 

This method is similar to method 5 except that ri* is for a blackbody and is based on 
instead of X^, where = X^/e , and the nonblackbody effect is accounted for by an apar- 
ent emissivity e . The apparent emissivity accounts for the multiple reflections among 
the fin and tube surfaces by considering the fin and tubes to be an isothermal cavity of 
specified local surface emittance e . The apparent emissivity is a function of e (held 
constant in this report), f/r^, and Tg/T^, while rj* is a function of N^, and 

Tg/T^. This method is developed in reference 5. 


Overall Radiator Efficiency 

The heat rejection rates calculated by both the one-dimensional and two-dimensional 
analyses are normalized by dividing by an ’’ideal" heat rejection rate to form a so- 
called radiator or heat rejection efficiency; thus, 


Q 

7 ] 


Q. 


id 


(19) 


The advantage of comparing all the analyses on the basis of ?j is that the resulting values 
always lie between 0 and 1 (except method 1 in one case) regardless of the size or temper- 
ature level of the radiator section. The normalizing factor was arbitrarily chosen as 

«id ='‘’(4 +--o) 

that is, the ideal heat rejection rate is based on the inside wall temperature and the pro- 
jected area of the tube. 

An additional advantage for the use of is that the efficiency for the one- 
dimensional analyses can be expressed in terms of nondimensional parameters, and 
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these same parameters can be carried over to the two-dimensional analysis. 

Two-dimensional . - To facilitate comparison of one- and two-dimensional methods, 
the two-dimensional equations are nondimensionalized in appendix D in such a way that 
the usual one-dimensional parameters appear. It is shown that the two-dimensional heat 
transfer depends on (which are one-dimensional parameters) and also on 

^o/^in 

The two-dimensional efficiency 77213 3^® ^2D^^id t»y the ratio of 

equations (12) to (20), and in terms of the dimensionless parameters of appendix D, is 
given by 


’’ 2 D ~ 


/ ^=(7t/ 2) yiX=R +L 

-^R d<p- j 

1 3R ° / Z9Y 

sin"\t/f)(jl>/r ) •/x=Vr2-L^ 


dX 




( 21 ) 


One-dimensional . - The one-dimensional efficiency, denoted by 77 . for method i, is 
defined as terms of dimensionless parameters, these ratios become as 

follows: 


”1 = 



r 


o 


'in 


+ 1 


”2 = 




”3 = 



( 22 ) 


(23) 


(24) 
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( 25 ) 


V4 = 



-^5 = 



^6 = 



(26) 


(27) 


where 


e 


o 



RESULTS AND DISCUSSION 

Results were obtained for several cases divided into three groups as indicated in 
table n. Group A consists of several radiators that were designed for various power 
system applications and are near minimum weight design. These radiators cover a wide 
range of actual dimensional variables. Group B is a set of parametric cases obtained by 
choosing three values of each dimensionless variable in order to study the effect on one- 
and two-dimensional agreement of varying each parameter independently of the others. 
Group C is a set of off-design Brayton cycle cases. The basic design case is a near- 
minimum weight radiator for a Brayton cycle space powerplant. This basic design case 
was altered by arbitrarily varying the inside tube temperature and the outside tube radius. 


Radiating Effectiveness 

The first results to be presented are the comparison of the one- dimensional radiator 
effectiveness with the two-dimensional shown in figures 5(a) to (f). The comparison for 
method 1 is shown in figure 5(a). In all cases the one-dimensional effectiveness is too 
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TABLE n. - INPUT DATA AND PERCENT DIFFERENCE RESULTS 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 


10 

11 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 
’'15 


Mate- 




Physical dimensions 




Dimensionless parameters 


Percent difference 

rial 

Emit- 

Temper- 

Inside 

Outside 

Fin half 

Fin half 

Thermal 

Equiva- 

Conductance 

Fin -tube 

Radius 

Fin aspect 

Sink 

®1 

^2 

®3 

®4 

=5 



tance, 

ature, 

tube 

tube 

thick- 

length, 

conduc- 

lent 

parameter 

profile 

ratio. 

ratio. 

temper- 








€ 

inside 

radius. 

radius, 

ness. 

i. 

tivity, 

sink 



ratio. 

V^ln 

t/f 

ature 









tube 

’'in' 

*• 0 ' 

t. 

in. 

K 

temper- 


^0 

*/^o 



ratio. 









wall, 

in. 

in. 

in. 


Btu 

ature. 















Tin. 





(hr)(ft)(°R) 

Ts« 















°R 






°R 














Design cases (A) 


A1 

0.9 

748 

0.500 

0.695 

0.0084 

4.17 

A1 



707 

.602 

.801 

.0129 

5.20 

A1 



838 

.095 

.293 

.0025 

1.58 

Be 



1149 

.187 

.355 

.0078 

2.13 

A1 



1149 

.250 

.511 

.0106 

3. 07 

SS 



1149 

.187 

.354 

.0158 

1.06 

A1 



607 

.062 

.213 

.0034 

2.56 

Be 



1656 

.312 

.813 

.0572 

1.63 

Be 



1700 

.312 

.885 

.0501 

1.74 

Cb 



1664 

.250 

.673 

.0402 

1.46 

A1 



1125 

.086 

. 320 

.0184 

2. 37 

Cb 



2210 

.187 

. 381 

.0264 

.763 

Be 



1670 

.187 

.441 

.0258 

1.10 

1 (b) 

(b) 

1 (b) 

1 <"> 

1 (b) 

1 (b) 

1 "?> 1 


536 

756 

976 

1196 

1400 

1550 

1700 

1850 

2000 

2130 


1.8 

.9 

1.5 


0.0770 7.21 


111 

110 

108 

75.5 
112 

10 

110 

54 

51.5 
34 

90 

38 

49.5 


400 

1.001 

1.00 

6.00 

1.39 

0.00201 

0.535 

8.6 

-4.5 

-4.5 

-7.7 

-2.0 

-3.8 



.864 

.864 

6. 50 

1.33 

.00248 

.566 

8.6 

-3.4 

-3.4 

-7.7 

-1.3 

-3.2 



.705 

.705 

5.40 

3.07 

.00156 

.477 

8.6 

-4.4 

-4.4 

-7.8 

-2.3 

-3.7 



1.50 

1.60 

6.00 

1.89 

.00367 

.348 

11.2 

-3.5 

-3.7 

-5.7 

-.3 

-2.2 



1.54 

1.54 

6.00 

2.04 

.00347 

.348 

11.1 

-3.6 

-3.8 

-5.7 

-.4 

-2.3 



1.39 

1.38 

3.00 

1.89 

.0149 

.348 

19.6 

-3.3 

-4.7 

.1 

1.5 

.9 



.503 

.502 

12.00 

3. 41 

.00133 

.659 

2.6 

-3.6 

-3.7 

-12.8 

-3.3 

-5.6 



.499 

.487 

2.00 

2.60 

.0352 

.242 

25.4 

.7 

-2.2 

3.5 

2.0 

2.4 

0 

.740 

.719 

2.08 

2.67 

.0288 

0 

26.0 

.4 

-3.0 

3.5 

2.3 

2.6 

0 

.924 

.894 ' 

2.17 

2.69 

.0275 

0 

26.7 

.5 

-3.3 

3.5 

2.7 

2.8 

700 

.623 

.622 

7.40 ' 

3.74 

. 00774 

.622 

9.8 

-.5 

-.7 

-7.0 

.6 

-1.4 

0 

.804 

.781 

2.00 

2.03 

.0346 

0 

26.4 

-.1 

-3.4 

3.9 

2.4 

2.7 

0 

.569 

.561 

2. 50 

2. 35 

.0234 

0 

20.8 

-.9 

-2.5 

1.1 

1.4 

1.5 


Parametric cases (B) 


(b) 


(b) 

1. 

0 

0.992 

6. 

0 

2. 

0 

0. 0250 

0.4 

15.7 

1.8 

0.9 

-2.6 

3.5 

1.6 




5 

.498 

6. 

0 







12.7 

1.1 

.6 

-4.7 

2.0 

.6 



2. 

0 

1.97 

6. 

0 







19.5 

2.5 

.7 

-.3 

5.0 

3.0 



1. 

0 

.953 

1. 

0 







33.8 

-2.5 

-8.4 

7.0 

-.1 

1.7 





.997 

15.0 







9.5 

3.0 

2.6 

-7.2 

3.2 

0 





.996 

6. 

0 

1. 

4 





13.7 

-1 

-.4 

-3.8 

2.2 

-.4 





.987 



3.0 





17.4 

3.4 

1.9 

-1.7 

4.5 

2.6 





.984 



2 

0 

.05 



20.0 

5.6 

3.7 

.1 1 

6.4 

4.5 





.994 





.0167 



13.9 

.3 

-.4 

-3.8 

2.3 

.4 





.995 





.025 

.8 

16.3 

2.0 

1.0 

-2.1 

3.9 

1.8 





.991 





.025 

0 

15.6 

1.8 

.9 

-2.7 

3.5 

1.6 


Off -design Brayton cases (C) 


400 


0. 

999 

0.994 

6. 

0 

2. 

4 

0. 0107 

0. 

452 

12.8 

-0.7 

-1.2 

-4.6 

1.4 

-0.5 



.998 

12.0 

1. 

2 





6.9 

-.8 

-.8 

-9.2 

.1 

-2.7 



,989 

4. 

0 

3.6 





16.5 

-1.4 

-2,5 

-2.1 

1.4 

.3 



.997 

8. 

0 

1. 

8 





10.4 

-.4 

-.7 

-6.4 

1.1 

-1.3 



.992 

4. 

8 

3 . 

0 





14.8 

-1.1 

-1.9 

-3.2 

1.5 

0 


223 

.223 

6. 

0 

2.4 



.746 

9.2 

-.7 

-.8 

-7.3 

-.6 

-1.3 

,625 

.623 







.529 

11.4 

-.6 

-1.0 

-5.6 

.8 

-.8 

1. 

34 

1.34 







.410 

13.8 

-.8 

-1.5 

-3.9 

1.7 

-.2 

2.47 

2.45 







.334 

16.2 

-1.2 

-2.3 

-2.4 

2.5 

.5 

3.97 

3.90 







.286 

18.5 

-1.4 

-3.3 

-1.2 

3.1 

1.2 

5. 

38 

5.25 







.258 

20.3 

-1.5 

-4.0 

-.5 

3.5 

1.7 

7. 

10 

6.88 







.235 

22.1 

-1.5 ' 

-4.7 

.1 

3.8 

2.2 

9. 

16 

8.79 







.216 

24.0 

-1.3 

-5.5 

.7 

4.1 

2.6 

11. 

57 

10. 99 







.200 

26.0 

-.9 

-6.2 

1.1 

4.4 

3.0 

|l3.97 

13.15 





1 


.188 

27.9 

-.4 

-6. 8 

1.4 

4.7 

3.4 


is not an independent parameter in the two-dimensional method but is in some of the one -dimensional methods. 
^’Physical dimensions not applicable. 
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large, all the points falling above the line of equality. This is because method 1 neglects 
the temperature drop in the tube and uses the actual tube area without allowing for any 
interchange with the result that the calculated tube heat rejection rate is too high. Thus 
the worst points on figure 5(a) are for cases having a relatively small value of (as 

can be seen in table D), which makes the tube heat rejection contribution more important. 

Figure 5(b) shows the comparison for method 2. Here the agreement is quite good 
even though the temperature drop in the tube is neglected. 

The comparison for method 3 is shown in figure 5(c). Here the agreement is also 
quite good but no better than for method 2. Method 3 is the same as method 2 except that 
the temperature drop in the tube is taken into account. Thus, for the range of cases 
studied, the tube wall temperature drop seems relatively unimportant. 

The comparison for method 4 is shown in figure 5(d). The agreement is fair but there 
is considerable scatter. This is probably due to the attempt to account for the radiant 
interchange by using the constant factor of 0. 85. 

Figure 5(e) shows the comparison for method 5. Here the agreement is quite good 
and there is very little scatter. It is no better, however, than methods 2 or 3 which are 
simpler. There is also quite good agreement and very little scatter in figure 5(f), which 
shows the comparison for method 6. Method 6 is not significantly better than methods 2 
or 3 but it is the most complicated of all the one-dimensional methods. 

To summarize, one -dimensional methods 2, 3, 5, and 6 show good agreement with 
two-dimensional; whereas method 1 shows very poor agreement and method 4 shows fair 
agreement. 


Percent Difference 


The agreement between one- and two-dimensional calculation can also be evaluated 
by the percent difference defined as 


Ej = 100 





(28) 


The percent difference for each method and for each case analyzed is given in table n. 
Here the agreement can be seen in relation to the input parameters of the various groups 
of cases. It should be noted that a positive percent difference means that the one- 
dimensional method is predicting too high a heat rejection rate, and a negative percent 
difference means too low a rate. 

The percent difference data is summarized in table HI in the form of average abso- 
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lute values of percent difference and 
the maximum absolute value of percent 
difference. As in figure 5, it can be 
seen here that methods 2, 3, 5, and 
6 are in good agreement with the two- 
dimensional method. Ilie best meth- 
ods, on the average, are 2 and 6 , with 
values of average absolute percent dif- 
ference of 1. 8 and 1. 9, respectively. 
Both methods have a value of 5 . 6 for 
the maximum absolute percent differ- 
ence, which is an indication of the per- 
formance of a method at its worst. 

Thus it appears that for the average 
radiator configurations likely to be 
encountered, method 2 will give as close agreement with the two-dimensional as method 6 . 
Method 2 has the advantage of being simpler than method 6 . 

The heat rejection formula for method 2 (eq. (14)) contains only prescribed quantities 
except for 77 ^^^ which is a function of and T^/T.^^. Curves for are readily 
available in the literature (e. g. , ref. 9). Furthermore, T /T. can often be assumed to 
be zero and only one curve against X.^^) need be used. On the other hand, the for- 
mula for method 6 (eq. (18)), contains three quantities not directly prescribed T^, rj*, 
and e . The outside tube wall temperature T^ is calculated by assuming one-dimensional 
radial heat conduction in a radiating tube. The overall blackbody effectiveness rj* is a 
function of N^, f and T^/T^; Tg/T^ is often assumed to be zero and 77 * is obtained 
by interpolation from plots of 77 * against and (ref. 5). The apparent emissivity 

e is a function of e, £/r , and T/T: here T^/T^ is assumed to be zero and e is 
obtained by interpolation from plots of e against and e (ref. 5). 

Because of its simplicity and accuracy, method 2 is recommended for central fin- 
tube radiator calculations except for cases with very large r^/rj (say, greater than 4. 0). 
For these latter cases, method 3, which is only slightly more complicated (the temper- 
ature drop in the tube must be calculated), can be used. 

Effect of Input Parameters on Percent Difference 

The cases of group B in table n were analyzed in order to determine the effect of the 
independent d i mensionless parameters on percent difference. Case B1 is the standard 
case and its parameters have typical values as obtained in current near-minimum-weight 


TABLE m. - SUMMARY OF DIFFERENCE DATA 


Group 

Number 

Method 


of cases 

1 

2 

3 

4 

5 

6 

Average absolute percent difference 

A 

13 

15.8 

2.3 

3.3 

5.4 

1.7 

2.7 

B 

11 

17.1 

2.2 

1.9 

3.3 

3. 3 

1.7 

C 

15 

16.7 

1.0 

2.9 

3.3 

2.3 

1.4 

Overall 

39 

16.5 

1.8 

2.8 

4.0 

2.4 

1.9 

(A, B, and C) 








Maximum absolute percent difference 

A 

13 

26.5 

4.5 

4.7 

12.8 

3.3 

5.6 

B 

11 

33.8 

5.6 

8.4 

7.0 

6.4 

4.5 

C 

15 

27.9 

1. 5 

6.8 

9.2 

4.7 

3.4 
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Independent dimensionless parameters 

Figure 6. - Effect of change in independent dimensionless parameters on percent difference. 


design studies at Lewis. The other cases were obtained by holding all parameters but one 
constant and taking high and low extreme values spanning the range of current interest. 

For example, case B1 has a of 1. 0, case B2 of 0. 5, and case B3 of 2. 0 with the 
other parameters being identical for all three cases. This results in three data points 
for each parameter. The results are shown in figure 6 for one -dimensional methods 2, 

3, 5, and 6. It should be pointed out that the level of the curves (i. e. , their location rel- 
ative to zero percent difference) is the result of the particular set of parameters chosen 
for the standard case. 

In looking at all the parts of figure 6, it can be seen that, in general, variation in 
IL/t^ has the greatest effect (particularly at small variation in t/H has the next 

greatest effect, and variation in 9^ has practically no effect. Another general observa- 
tion is that all four methods exhibit essentially the same trends except for the i/r^ curve 
for method 6. 


CONCLUSIONS 

Results of the two-dimensional analysis of the heat transfer in the cross section of a 
central fin-tube radiator have been used to evaluate several one- dimensional methods of 
calculating the heat rejection rate of such radiators. Most of the one-dimensional methods 
give good agreement with the two-dimensional; however, they are of varying degrees of 
complexity. In view of the good agreement of the method (number 2) that neglects tube 
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wall temperature drop and accounts for radiant interchange simply by using the projected 
area of the tube, it seems unwarranted to use more complicated methods. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 30, 1966, 

120-27-04-36-22. 
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APPENDIX A 


SYMBOLS 


A area of radiator surface per unit 
axial length 

B radiosity, sum of emitted plus re- 

flected energy leaving surface 
per unit time and area 

M 

dA* particular element of area 

Ej percent difference of method i 

Fj _2 view factor, fraction of radiant 

energy leaving surface 1 that 
strikes surface 2 

H total energy incident on a surface 

per unit time and unit area 

k thermal conductivity, summation 

index (appendix C) 

L £t/£^ = t/£ 

£ fin half length, see fig. 2 

N„ blackbody conductance parameter, 

C Q O 

Q calculated heat rejection rate 

Qid ideal heat rejection rate, eq. (20) 

q energy radiated per unit time and 

per unit area 

R rt/£^ 

r radial coordinate in tube 

S arbitrary surface representing 

surroundings 

T temperature 


Tg equivalent sink temperature 

t fin half thickness, see fig. 2 

X xt/£^ 

X horizontal coordinate in fin 

Y yt/£^ 

y vertical coordinate in fin 

0! absorptance 

€ emittance 

e apparent emissivity 

Tj overall efficiency 

T?in fin efficiency based on 

r\^ fin efficiency based on 

7j* overall effectiveness, ref. 5 

Tjrp overall effectiveness, refs. 4 and 6 

9 dimensionless temperature, T/T.^^ 

X conductance parameter, eoT|^£ /kt 

p reflectance 

CT Stefan-Boltzmann constant 

(p angular coordinate in tube, see 

fig. 3 

o) overrelaxation parameter, appen- 

dix C 

Subscripts: 
b fin base 

2D based on two-dimensional analysis 

i based on one- dimensional method i, 

i = 1, 2 ... 6 
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u 


i, j,k indices, appendix C 

in inside tube wall 

L lower limit of visibility 

o outside tube wall 

S referring to area of arbitrary sur- 

face S 


1 

2 


upper limit of visibility 

tube 1 
tube 2 


Superscript: 

m iteration number, appendix C 



APPENDIX B 


DERIVATION OF EQUATION FOR INCIDENT RADIATION 

The total radiation incident on an element dA* of the fin or tube is made up 

of radiation from other parts of the fin or tube and radiation from the surroundings, hi 
equation form, 

**dA* " ^A~dA* ‘^S-dA* 

where is the heat per unit time leaving all of the fin and tube surface A and 

striking a unit area of dA*, and is the heat per unit time leaving the surround- 

ings or sink which may be represented by an arbitrary surface S (fig. 7) and striking a 
unit area of dA*. The q’s are given by 


‘iA-dA* ■ 


^A-dA* 

dA* 


and 


_ ^S-dA* 
<lS-dA* = 


(B2) 


(B3) 


where Q^^_(Ja* heat per unit time leaving all of A and striking all of dA*, and 

Qg_dA* is that leaving all of S and striking all of dA*. The Q's are given by 


*^A-dA* ■ 



(B4) 


and 


^S-dA ■ ‘^^S-dA* 


(B5) 


where is the radiosity, that is, the total energy emitted plus 

per unit time, is the view factor from dA to dA*, T^ 

temperature of the surroundings, S is the area of surface S, and 


reflected leaving dA 
is the equivalent sink 
dFg_dA* is the view 
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factor from S to dA*. From the reciprocity theorem for view factors, 


^ ^^dA-dA* “ *^^dA*-dA 


(B6) 


and 


® ^^S-dA* ■ ‘^*^dA*-S 


(B7) 



Furthermore, all the energy leaving dA* 
strikes either part of S or part of A (see 
fig. 7) so that 


^dA*-S = 1 



(B8) 


Combining equations (B6), (B4), and (B2) yields 


/ 


^A-dA*" / ®dA^^dA*-dA 


(B9) 


Combining equations (B8), (B7), (B5), and (B3), yields 


1s-dA* = 


{■/ 


^^dA*-dA, 


(BIO) 


Putting equations (B9) and (BIO) into equation (Bl) yields the expression for the total en- 
ergy incident on dA* per unit time: 
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APPENDIX C 


NUMERICAL METHOD OF SOLUTION 

The nonlinearity of the boundary conditions and the geometry preclude the possibility 
of a closed form solution of the problem outlined in the Equations and Boundary Conditions 
section (p. 4). Hence, it was necessary to attack the problem numerically. The numer- 
ical solution was obtained by replacing equations (2), (3), and (10) by their finite- 
difference analogues at points of a grid as shown in figure 3 (p. 4). The solution of the 
partial differential equations satisfies these difference equations except for an error term 
which vanishes when the grid spacing approaches zero. This in turn is reflected in an 
error in the temperatures at the grid points which likewise vanishes as the grid is refined. 

The set of points at which the temperatures are obtained consists of the points of in- 
tersection of the polar grid lines in the tube and the points of intersection of the rectangu- 
lar grid lines in the fin (fig. 3). The spacings of the rays and circular arcs in the tube 
and verticals in the fin are prescribed. The arcs in the tube and the verticals in the fin 
are uniformly spaced. For the spacing of rays in the tube, however, the tube is divided 
into two regions of uniform ray spacings; the ray spacings near the fin are finer than the 
spacings away from the fin (as illustrated in fig. 3) since the temperature gradients are 
more severe near the fin. The spacing of the horizontals in the fin is determined by the 
ends of the rays in the tube at the intersection of the tubes and fin and is thus not uniform. 
After this main grid is specified, a dual grid is constructed by bisecting the lines of the 
grid just described. This results in a closed cell about each point of the original grid. 

The boundaries of the cells are the dashed lines shown in figure 8. 

The finite- difference equations for the temperatures are obtained by integrating 
Laplace's equation over cells of the dual grid, replacing the normal derivatives arising 
from the use of the divergence theorem by difference quotients exactly as on pages 17 and 
18 of reference 11. In the present case, however, the cells of the dual grid are not all 

rectangles as they are in reference 11, but some are 
annular sectors in the tube and some are irregularly 
shaped cells in the transition region, one of which is 
shaded in figure 8. The result is that the difference 
equations in the tube have coefficients modified by 
lengths of arcs and chords and those in the transition 
region entail six points instead of five. 

The finite-difference equations for the radiosities 
oi ?n"and ^re obtained by substituting a simple quadrature for- 

mula for the integral in equation (10), so that it reads 
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where the summation is carried out over those elements dAj^ that are visible from ele- 
ment dAj. 

The result of the discretization is two sets of difference equations: one set is derived 
from equations (2) and (3) at interior points and from equations (2) to (7b) at boundary 
points, and a second set is derived from equation (10). The first set is a system of equa- 
tions for the temperature interior to and on the surface of the tube and fin with the values 
of radiosities as input to the equations at the surface points. These latter equations were 

linearized by replacing T^ by +4 [t^*^”^^] [t - wherever it appears. 

Here the superscripts refer to iteration number, and m = 1 corresponds to the initial 
approximation. The second set, typified by equation (Cl) is a system of linear equations 
for the radiosities with the values of surface temperatures as input. Initially, very 
coarse approximations were made to both temperature and radiosity. Then one iteration 
of the block successive over relaxation method was carried out on the first set of equa- 
tions. This single iteration for the temperatures was followed by as many iterations as 
needed for convergence of the Gauss-Seidel method on the linear system of equations for 
the radiosities (very few iterations were needed for this step). This process was re- 
peated, alternating between the two sets of equations until convergence was obtained. It 
was found that comparison of the net heat radiated by the fin and tube surfaces with that 
flowing into the inside surface of the tube supplied a practical criterion for the conver- 
gence of the iterative process. 

The block successive over relaxation method mentioned previously was like that 
described on pages 19 and 20 of reference 11, at least in the fin. In the tube, the itera- 
tion proceeded along similar lines except that here temperatures along a given ray were 
solved for simultaneously. The iteration started with the vertical ray (at (p = tt/ 2), pro- 
ceeded to each ray in turn until the horizontal ray {(p = 0) was reached. Then the itera- 
tion proceeded to the verticals from left to right in the fin. 

In the iterative method mentioned previously, there is an overrelaxation parameter 
CO, which is free to be chosen subject to the restriction 0 < co < 2. It has been shown 
(ref. 12, ch. 4) for systems of linear equations that there is an optimum value of u> for 
which convergence is fastest. Due to the nonlinearity of the problem involved herein, no 
rigorous analysis exists for the determination of a corresponding optimum value for co. 
Numerical experiments carried out at the Lewis Research Center, however, indicated 
the validity of the concept of an optimum value, and values of co obtained in these experi- 
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ments were used for the production runs. A typical result is o) = 1. 93 for 690 points 
leading to convergence of the relative error in the heat balance to less than 1. 0 percent 
in 180 iterations. 



APPENDIX D 


NONDIMENSIONAL PARAMETERS AND EQUATIONS 


In order to facilitate comparison of the two-dimensional results with one-dimensional 
results, several dimensionless quantities are formed. This set of parameters results 
from forcing the conventional one- dimensional conductance parameter to appear; 


e 




e 


s 



X 



»ln- 


kt 



Y 






(Dl) 


Putting these parameters into the dimensional equations and boundary conditions results 
in the following: 


In the tube: 


lii = 0 

^ 3R^ d(p^ 


In the fin: 


= 0 

3X^ 3Y^ 

2 

Along the inside tube wall where R = Rj^^^ = = (rjj^/£)(t/£): 


(D3) 
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= 1.0 


(D4) 


Along (p = 0 and <p = v/2 in the tube: 

/de\ 

y^)(p=0 

Along Y = 0 and X = (£ + = (1 + r^i)t/S.\ 


= 

W/(p=(7r/2) 




= 0 


X=(l+rQ/£)tA 


On the surface of the tube where R = = r^/S. = (r^/£)(t/£): 


ra^' 

^9R> 


= Xin(0 - ^ 


R=R. 


2 2 

On the surface of the fin where Y = t /£ : 


— J = ^ 

WY.t2/,2 


(D5) 


(D6) 


(D7a) 


(D7b) 


where 



Thus, it can be seen that the dimensionless temperature 6 depends on the parameters 

X. , r. A, r /£, t/£, 6^, and e. In this report, e will always be fixed at 0. 9 and can 
in^ in' ^ o' ^ ^ S' 
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be eliminated from the list of variables. Of the remaining parameters, r^/S. (often 

written £/r ), and 6„ are the commonly used one-dimensional parameters and r. /£ 
o s in 

and t/£ are two-dimensional parameters, bistead of it is more informative to 

use which is the ratio of r^£ to To summarize, in this report the 

parameters determining the two-dimensional heat transfer in a fin-tube radiator are X^^, 
r^/S. (or £/r^), 0^, r^r.^^, and t/£. 
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